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文章 编号 :1005-3085(2011)02-0165-11 
一 类 基于 综合 害虫 管理 策略 的 捕食 系统 的 动力 学 性 质 * 


王 玲 书 22， 徐 瑞 2?， 汉 光辉 ? 
(1- 河北 经 贸 大 学 数学 与 统计 学 学 院 ， 石 家 庄 050061; 
2- 军械 工程 学 院 应 用 数学 研究 所 ， 石 家 庄 050003) 
摘 要 : 本 文 基于 综合 害虫 管理 策略 ， 研 究 了 一 类 具有 脉冲 效应 的 捕食 模型 并 对 其 动力 学 性 质 进 行 了 分 
Pi. 根据 Floquet 乘 子 理论 和 脉冲 微分 方程 的 比较 定理 ， 证 明了 当 脉 冲 周 期 小 于 某 个 临界 值 时 ， 
系统 存在 一 个 全 局 渐 近 稳定 的 害虫 根除 周期 解 ， 给 出 了 系统 持续 生存 的 充分 条 件 . 数值 分 析 进 一 
步 说 明了 系统 的 动力 学 性 质 及 综合 害虫 管理 策略 的 有 效 性 . 
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1 引言 


农药 的 长 期 使 用 ， 使 害虫 产生 了 抗 药性 ， 再 加 上 害虫 的 高 出 生 率 及 气候 等 原因 ， 虫 害 很 容 
易 爆 发 . 例如 ，2000 年 美国 密 执 安 州 大 豆 蚜 虫 的 泛滥 ，2002 年 中 国 -一 些 地 区 的 蝗灾 等 都 造成 
了 严重 的 经 济 损失 . 因此 ， 害 虫 控 制 是 关系 到 经 济 发 展 的 一 个 非常 重要 的 问题 . 近 二 十 年 来 ， 
出 现 了 许多 控制 害虫 的 方法 ， 其 中 ， 综 合 害虫 管理 上 3] 是 使 经 济 损失 达到 最 小 的 一 种 非常 有 效 
的 方法 ， 它 主要 包括 生物 控制 和 化 学 控制 两 种 方法 .生物 控制 中 通过 释放 天 敌 (捕食 者 ) 或 利用 
病理 来 控制 害虫 ( 食 饵 )， 是 近年 来 备 受 重视 的 生物 防治 的 一 个 领域 ， 因 其 可 以 避免 化 学 药剂 带 
来 的 问题 而 且 益 受到 重视 ， 化 学 控制 是 通过 喷洒 杀 虫 剂 来 控制 害虫 ， 它 能 使 害虫 数量 迅速 减 
少 ， 尤 其 当 害 虫 数 量 太 大 释放 天 敌 不 足以 控制 害虫 或 考虑 到 释放 天 敌 的 成 本 时 ， 必 须 使 用 杀 虫 
剂 来 控制 害虫 .实践 证 明 ， 综 合 害虫 管理 比 任何 一 种 传统 的 方法 都 更 有 效 . 

在 经 典 的 捕食 模型 中 ， 人 们 总 假定 每 个 捕食 者 具有 相同 的 捕食 和 生产 能 力 . 这 个 假设 对 许 
多 动物 来 说 似乎 是 不 切实 际 的 . 在 自然 界 中 ， 物 种 的 增长 常常 有 一 个 成 长 发 育 的 过 程 ， 如 从 幼 
年 到 成 年 等 ， 而 且 在 其 成 长 的 每 一 个 阶段 都 会 表现 出 不 同 的 特征 . 因此 ， 考 虑 具有 阶段 结构 的 
种 群 模型 更 具有 实际 意义 .近年 来 ， 具 有 阶段 结构 的 种 群 动力 学 模型 已 引起 了 许多 学 者 的 关 
注 . 在 文献 [4,5] 中 ， 作 者 考虑 了 下 面 的 捕食 模型 

i(t) = rz(t) 一 az2( — 42820 


1+mz(t) ， 
h(t) = PG — (rı + d)y (0). (1) 
y(t) = dyı(t) — r2y2(t), 
其 中 z(t) 表示 食 饵 (害虫 ) ERA t REE, y(t), yo(t) 分 别 表示 捕食 者 (RH) 幼体 和 成 体 
在 时 刻 t 的 种 群 密度 ，o, a1, az, d, r, ri 和 2 为 正常 数 ， 其 具体 的 生物 意义 可 以 参考 文献 团 ， 
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对 于 系统 (1)， 由 文献 [4] 可 知 ， 平 衡 点 1(0,0,0) 是 不 稳定 的 ， 当 ra(r + d) > 如 和 时， 边界 


+mr 
平衡 点 BE2(5,0,0) 是 全 局 渐 近 稳定 的 ， 当 72(71 td) < Ait, RAER EER. 因此 ， 用 


一 次 性 引进 天 敌 的 方法 控制 害虫 并 不 有 效 . 本 文通 过 周期 投放 天 敌 和 喷洒 杀 虫 剂 对 系统 (1) 3E 
行 改进 ， 考 虑 如 下 脉冲 微分 方程 








i(t) = re(t) — az?(t) Mie 


y(t) = ezO) — (ry -d(t) ptz (nl—1)r, tznr, 
yalt) = dw (t) — raus (t) 


Ax(t) = —pz(t) 


Ay (t) 2 -piy (t) e t9 (n L- 1r, (2) 
Ay2(t) = —p2y2(t) 

Az(t) = 

Ay (t) =u, (tmr, 

Ayx(t) = 


其 中 和信 z(t) = z(t*) — r(t), Ayit) = w(t*) yilt) ræ, ne 2Z+ = {1,2}, 0< 
p, pi < 1 分 别 表 示 在 1 = (n+l 1)r EtA t F2 k HC 6 E B AR E T S EG 
Bl, u > 0, i=1, 2 EFt — mr 时 刻 释放 天 敌 的 数量 . 


2 ”害虫 根除 周期 解 的 全 局 渐 近 稳定 性 


S f= (fy, fo, fs) 是 系统 (2) 前 三 个 方程 的 右 端 映射 ， 显然 了 的 光滑 性 保证 了 解 的 存在 唯一 
性 . 设 z(t) = (4), yi (t), y(t) 是 系统 (2) 的 解 ， 容 易 证 明 下 面 引 理 成 立 . 


引 理 2.1 系统 (2) 的 满足 初始 条 件 z(0+) > 0 的 解 ， 于 是 对 所 有 的 t+ > 0 都 


引 理 2.2 Sp 的 所 有 正解 是 一 致 最 终 有 界 的 ， 即 存在 一 个 正常 数 M， 当 t 充 分 大 时 ， 
有 z(t) € M, yilt) < M (i = 1,2). 


考虑 (2) 的 子 系统 


y(t) = —(rı + dt(t) 
yo(t) = dyi (t) ~ r2y2(t) 
y (tt) = (1 - nuit) lisasin (3) 
yo(t*) = (1— p2)y2(t) 
yi(tt) = u, + Y(t) | 

t = 2T. 


je 上 天 mT， 


yo(t*) = u2 + ya(t) 
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引 理 2.3 系统 (3) 存在 一 个 全 局 渐 近 稳定 的 正 周期 解 (y*(t), yt) HP 


uie (rid) n7) 


cujye MITEN Lm 
1—(1-pije- C137 nr «t& (n t D7, 


1-— —(Gi-td)(t-nr) 
pU ena (nt)r«t£(n4 Ur 














| dup. L E e (r1 t97 +(p2-pi)e C1td-r2)lr-r27)] —rz(t-n7) 
1-(17pa)e "27 (ri*d-r3)üi-ü-p)e (1*27)ü-ü-pa)je-737) |? 
dujie-C1t*4)(t-n7) 
7 ira (perma "T< tS (n+D), 
*(t)— v2 (1—p2) dui(1—p2)(p2—p1)e (71 *4—772)!7 
v(e) [e + (rı+d-r2)(1—(1—p1)e 701 *07)(1—(1—pa)e7727) 
+ dui(1l—pz) e-r2(t-nT) 
(ri 十 dg 一 r2)(1 一 (1 一 p2)e-727) |È 
dui(1—pi)e ^ (r3 d)(t—-n7T) 
atdi nA- doe Uie (n - Dr «t € (n4 1)r. 
证 明 在 区 间 nr «t € (n 十 1)+ 上 求解 (3)， 可 以 得 到 
wi (t) = w(nr*)e 0110070 
y2(t) = ya(nr*)e 207^? nr «t € (n4 )r, 
$ dn(n77) [e r2(t-nT) .— e (r14-d)(t r) 
y(t) = y (n Dr* Je Ciroe-emr) 
(5) 
y2(t) = ya((n + D)r*)e 726-007) (n 4- Dr « t € (n * 1)r, 


十 dnd [e 一 r2(t 一 (mn 十 DT7) _ errorum) 
y(nr*) = wi t yi(nr), 
y((n-D)r*) 2 (1 p)w((n-0)7), i-1,2. 
4 y(nr*) = y((n-- 1)*), i9 12， 代 入 上 式 可 以 得 到 系统 (3) 的 + 周期 解 UO, yO). 
进一步 ， 由 (5) 可 以 得 到 人 迭代 映射 
yi(nr*) = ui + (1 — pii ((n — 1)rt)e 2r, 


y (n Dr*) = (1 — pi)y (nrtje t9, 








€ d(1— —T2T a-(ritd—-r2)r—rT97T 
ys(nr*) = uz + (1 — pa)e7"?7 ys ((n — 1)r*) + Zee EL ————w((n- 1r*) 
d(e- 720-07. e-(144)0-07 
十 E 3yi ((n 4-1 — 1)7*), 
—ralr o7 (Gi d)ir 
yi ((n  Dr*) = (1 — pz)us(nr*)e7ritn + Hee 777-9 i (nr), 


由 此 可 以 得 到 


ui 1 = (1 — p1)^e-n(ri-td)r 


| * ne—n(ritad)r 
1-(ü-p)e mide + yi(07)(L— p)" e- "rinm. 





yi (nr*) = 


wll ~ (1 -pi)"e etr) 
1- (1 zm pile^ (rı+d)r 





y((n+1)r*) = (1— pi)e Siom +y (0*)(1 二 
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d(1 ~—p2)(e 27 —e (ri--d—r2)lr nr) 


+) 
Ti 十 Q 一 72 





y2(nTt) = [uz + Vi(nT 





e-72(1-07 — e-(ritad)(1-D)rj1— (1— paye7"r27 
T1 - d — r3 | 1 — (1 — p2)e™"27 


+ ya(0*)(1— pa)" e727, 


d(1 — pa)y (nr*)(e7?" — e) 
Tj 十 Q 一 T2 


+ dyi ((n+i)r* 





y2((n + D)r*) = +1- pa)yo(nr* je", 


这 样 就 可 以 得 到 
lim, (n), wo (0) = (vi (t). u3(0), 


即 系统 (3) 的 周期 解 是 全 局 渐 近 稳定 的 . 
定理 2.1 设 z(t) = (z(t) y (t), yo (t)) 是 系统 (2) 的 具有 正 初 值 的 解 ， 假 设 


aai v 1 
T >(t)dt < 1 
2* um ve TES 


成 立 ， 则 (0, y(t) yt) 是 全 局 渐 近 稳定 的 ， 其 中 


* — 1 dui —ralr _ 2 —r2T 
] i6 7 za a emn mre oem) 


2n pa)(pa pi) e7?7)e-( | 
(ri -d— T3)(1 - (1 — pi)e-(*97) 

















du, (1 22 pie 1t = (1 a pije (0*97) 
(rı - d)(ri +d — r3)(1 — (1 — p1)e-*27). 


证 明 首先 证 明 局 部 稳定 性 . 4 x(t) = elt), zat) = wi(t) yit), r(t) = yt) - v (t): 
则 系统 (2) 的 线性 化 系统 为 





i(t) = (r ~ asy (t))oi (t) 
g(t) = azy z(t )z1(t) 一 (r1 + d)x»(t) tz (n +l- 1)r, txnr, 


da(t) = dz2(t) — raxa(t) 

zi((n +l- 1)r*) = (1 7 pymi((n +1- 1)7), 
z5((n --1— 1)r*) = (1 — pi) ((n 十 — 1)7), 
z3((n - 0— 1)r*) = (1— pz)zs((n +l — 1)7), 
zi(nT+) = zi(n7), i— 1,2,3, 

定义 
r — a ys (t) 0 0 1-p 0 0 
A(t) — azt)  —-(n4d 0 |; B= 0 1-p 0 


0 d 一 72 0 0 l—po 
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根据 脉冲 微分 方程 的 Floquet 乘 子 理论 4 ， 如 果 单 值 矩 阵 


(1 — p)elc (r—aiy; (t))at 0 0 
M(r) = Belo A(t — * (1 — p)e- 0 *07 0 
* * (1 — pao)e-?7 


的 特征 值 的 模 小 于 1， 则 周期 解 (0, yr (0), y2 (0) 是 局 部 渐 近 稳定 的 .实际 上 ，M(7) 的 特征 值 
0< àz —(1— pe *?7 <1, 0<XM=(1-p2)e ??* <1, 


因此 当 
0< 和 1 一 (1 — pelo (r-aiy; (t))dt «1 


Wl. (0, y (t), y (0) 是 局 部 渐 近 稳定 的 . 
其 次 ， 证 明 (0,yf(£5),y3(0) 是 全 局 吸引 的 . 由 式 (6)， 可 以 取 <s > 0， 使 得 


c Ê (1 — pjel/é (0-3 (0-9)8t < 1, 


由 引 理 2.3 和 脉冲 微分 方程 的 比较 定理 有 ， 当 t 充 分 大 时 ， 不 等 式 yo(t) > y5(t) - eor. 不 妨 
假设 对 所 有 的 t > 0 上 式 成 立 ， 这 样 由 (2) 可 知 


| ilt) < z(t) fr — (y(t) —e)], tz n+l Dr i 


z(t*) = (1— p)a(t), t=(n+l-1)r, 
E ((n+l-1)r, (n+ D] 上 求解 上 式 ， 可 以 得 到 
z((n-- Dr) < z((n-1— 1)r*)efechimo 0-192 0-9) 

= z((n +l- 1)r)(1 — p)ele (7*2 (079t 

= z((n-4-1— 1)7)o, 
由 此 可 以 得 到 z((n - D)r) < z(ir)c"， 于 是 有 

„im z((n DT) = 0. 

另 一 方面 ， 当 ie ((n--L— Dr, (n - Dr]Hf, 0 < s(t) € z((n-1—1)7)0 — pe", 所 以 
当 t +œ}, z(t) 一 0， 由 此 可 知 存在 6， 使 得 t > To Ef, 0 < z(t) < e. 


不 妨 设 t > 0 时 , 0 < z(t) < e， 于 是 由 引 理 2.2 和 (2) 可 以 得 到 g(t) < agMe — (rı + 
d)yn(t). 考虑 


gi (£) = az2Me — (rı + d)y (t) 
ya(t) = dy (t) — r2y2(t) 
yi((n 十 1 一 1)r*) -(1 — pi)yi((n 十 1 一 1)r), 


yi(nr*) = ut y(nr), 1-12, 


je tnr, 
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对 任意 初 值 (yi(0+),ya(0+))， 求 解 (8) 可 以 得 到 ， 当 nr <t< (十 Dr 时 ， 有 





y(t) = yt) + (1 — p1)” n (0*) - i e 0 


pi(e1*2n7 (1p, )*)e- C16 707) | 
, 





Tid 


十 £22 M 1 a (1 " p1)"e- (1*2 i-ü-pie-(1*27 


li 


yt(t) 十 fient d foe—"2(t—n7) + ca3 Mdfs — dy1(0+)(1—p1)" e-(ritd)t 


y2 (t) rı+d-r2 


M (ncDr «tx (no l)r it, # 


y(t) = yi + 0-0)! [ui (0*) - Tuc et 


pie" (e(r1 *)n7 L(1-p,)nt!e- Cri t27)g- (ri tt 
1-(1-pi)e- 1*7 ; 


dyi(0*)1 — p)" ^ (ea 
Ti 4-d— r2 , 








ete [i-a carte 


yalt) = y3(t) + gie "2t + gae 7"? + eag Mdgs 





其 中 
— En" n m 
fi = (1 - p2) [us (0*) I- 35): 


duje-"?7(1 i p2)(1 ER gceerd- noma — (1 S CUR TE 
(r1 - d — r3)(1 — (1 — p1)e-1*297)(1 — (1 — pa)e77?7) ' 




















f E dyi (0*)(1 —p)*e Grtdwr f x e Direrd)r- mr 23 p2)(1 " (1 a paj"e-r2n7] 
Am Tı 十 Q 一 72 1— (1 — p2)e "27 
eh eritd—ra)lr | e-r2a(1-Dr | 
1 e-(ri*d)(t-nr7) e 72(t-n7) 
= + 
? rə(rı +d) (rı+d)(rı+d-rz2)  ra(ri-- d — r2) 
-ra(t-nT) 一 1- me 一 72t —ra(1-0)r —(r1i--d—ra3)ir—r2T 
mi e ( pa)"e a IE " e 
1— (1 — p2)e7"? rə(rı +d) (rı +d)(rı +d- r2) 


e 一 r27 d e-72(t-n7) E (1 ET p2)^e-"t 
ra(rı +d — 可 | (r1 *- d)(r1 +d — r2) (1 — (1 — p1)e-(1*27)(1 ~ (1— p3)e77?7) 


x [a „2 (1 _ pi)” e (^-D(-dr)( T pije (1*0 -07 ep pi)e 1*9] 


x (a E p3)e "?* (1 x ec Cxbaera)tm) 再 (1 o pije Uit (g- r8 0-07 E ea+d(G -07) 





Er (1 m e- (1 *dtry(e-r2G Dr e (rı+d)(1 07))， 


n u 
à - (- 5" n0) d er 


du; (1 Él pı)(e77277(rı+d-r2)lr TE e-(rutdinrya s (1 — pi)? 1e- 0+8)7) 
(rı +d — r2)(1 (Le pe terere (1 — p2)e tan) | 
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ides dyi(0*)(1 — p3)(1 — 91)" 
T1 十 Q@ 一 72 








e—"27(1 一 ea+e e) +e ro(1—l)r =e (ri+d)(1 | 


1 
| * 1 — (1 — pz)e-7727 





1 e- (rid) G7) e 727 (Q-H07) 


L (1 p erri onn) 
gs = (1+e Morum F (rı +d)(rı+d—r2) r2(rı+d-r2) 


+ (1— p3)* 








一 7r2(t 一 (PDT7) _ g- rid) (t- (07) e 720-07 e^ (ri d-r2)lr—raT 
十 
(rı 4- d)(r1 - d — r2) rs +d) (rı+d)(rı+d-r2) 
e 7?7 (1 e-(li+alir)(e-r20-07 E eri+dad-07) 


i T2(ri c d- r3) 二 (rı - d)(r1 +d — r2) 


ow acOibd)ü-Dr Aa —(ri4-d)r 
+ (1—pie (1 ~ pi)e ) 





(1 (1 E pi) le (n-0mt*dryie r2(1—l)r _ e (rı +d)(1—l)r + (1 Lm p2)e "7?" (1 I. e (itd rar 
A ` 
(rı +d)(rı +d — r2)(1 — (1 — p1)e71*97) 





RE EROAK, WAAR y(t) < yilt) +e vlt) < yzl) +e. 
男 一 方面 ， 由 (2) 可 以 得 到 加 (t) > 一 (ri 十 d)y(t). 于 是 由 引 理 2.3 和 比较 定理 可 知 ， 


当 t 充 分 大 时 ， 有 h(t) > (一 81, wt) » (t ^ &. 
由 以 上 分 析 可 知 


Jim. (y(t), v2(t)) = (vi (©), (0)). 


3 ”系统 的 持续 生存 


定理 3.1 如 果 





TT— a f yz(t)dt > ln (9) 
0 


1 
l-p 
成 立 ， 则 系统 (2) 是 持续 生存 的 . 

证 明 由 引 理 2.2 可 知 存在 常数 M > 0， 当 t 充 分 大 时 ,z(t) € M, y(t) € M (i = ,2), 
不 妨 设 t > 0 时 ， z(t) < M, w(t) < M. 因此 ， 要 证 系统 (2) 是 持续 生存 的 ， 仅 需 证 明 存 
在 常数 m; > 0( = 1,2,3), Ert) > m, y(t) > ms, y(t) > ms. B (2) T gt) > 
一 (71 十 d)yi(t)， 则 由 引 理 2.3 可 知 ， 对 任意 的 e > 0， 当 t 充 分 大 时 ， yilt) > yilt) ^e. 因此 ， 
当 t 充 分 大 时 ， 有 


ui(1 — pi)e (1*907 ^ 
y(t) 2 1- - pje ar — E = m». 
类 似 可 得 ， 当 上 充分 大 时 ， 有 
y(t) > u»(1 — p2)e "7 m ^ Re 





= 1-(1- pz)e7"27 


下 面 只 需要 证 明 存 在 一 个 正常 数 m1， 当 t 充 分 大 时 ，z(t) > mi. 我 们 分 两 步 来 证 明 . 
步骤 1 由 式 (9)， 取 


rr — a fj yi (t)dt — In 


0 < 
Sn (a + ayaogdMC}r 
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和 充分 小 的 sz， 使 得 


n 全 (1 e Djer7-m4(o+alo2dMC)T -a1 fo (và (£)2-e2)dt > 1, 


z 1 | 1 
| 1- (1 — pa)e77?7 Lr2(ri + d) 





C 


" e- (i t)r a de?" ((1 _ p2)(1 — e+d-r2)i7) E (1 — pie City (er20 pr 2 p PED -DTY 
(rı 4- d)(r1 +d — r2) 
e^7727 | 


z r2(ri +d— r2) 


下 面 证 明 一 定 存在 一 点 五 > 0, 1 x(t) > ma. T 
últ) < aym4M — (rı + d)yı (t). 


考虑 
0) = azm4M — (rı + di (t) 
Ya(t) = dyi (t) — r2y2(t) 
yi((n+1-1)r*) = (1 -piy ((n +1- 1)7), 
yi(nr*) = u, +yilnT), i=1,2. 


KAF (8) 可 以 得 到 ， 当 t 充 分 大 时 ， 下 面 不 等 式 成 立 yo(t) € y5(t) + maa2dMC + e2. 
设 te ((n 十 DD)7, (nn 十 1 十 1)7]， 由 此 可 以 得 到 


上 es tnr, 


| i(t) > z(t)|r — am4 — ai(y3(t) + maa2dMC + e). tz (n4), 
z(t*) = (1 — p)z(t), t= (n-- Dr. 


Hn +r >to EKE ((n +r, (n - L4 1)r] ERRERA 


z((n+l+1)r) > z((n+l)r)( — pelo -amaa (y (t)+maazdMC+ea)]dt 


= z((n +1)r)n, 


FEMk>oolt, z(n-- 1-4 k)r 2 z((n - D7)n* 一 00， 这 与 z(t) 的 有 界 性 矛盾 .因此 ， 存 在 
一 点 五 > 0， 使 得 z(t1) > m. 
步骤 2 ”如 果 上 之 五 时 ， 均 有 z( > ma， 则 结论 成 立 ， 否则 ， 令 


t = inf {z(t) < ma}, 
妇 分 为 脉冲 点 和 非 脉冲 点 两 种 情况 : 
1) AKIPA. QE = (ni 十 DT, ni € Zo Walt) z m, t e [tt] H(i- pm < 


z(t**) = (1 一 p)z(#*) < ma. 选取 ma, na € 2Z4+， 使 得 


1 E 
NT > In " 
T2 C1 





[15 secs PASA, 
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其 中 
01 — T —am4 —a4M « 0, 
U2 dui 1 
Cı = (1 — p2)™ 0 
1 ( p2) (vt ) 1— (1 二 pa)e-r27 ri -d—r531— (1 | 
d(1— p1)^ ^! (1 — p2) ui 
十 (ri gg rj) m (1 — p2)e—"27) (v. (o* )- 1- (1 Ense 


x ((1 — pae 27 十 (po — p1)e (1 td-7r2)tr—rar z (1 - pie i97), 

4 T = (no - na)r, W—5ETEfEt, € [t,1* - T], E elt) > m, WW Hte [i ea 

TIM, iji) < amaM — (red) (0). 类 似 于 (8) 可 以 得 到 ， 当 te [P emori T] M, 有 
|y2(t) 到 y»(t) 一 ma4a2dMC| ze Patar. 


MT yalt) € y3(t) + mad MC + css， 同步 骤 1 类 似 ， "If z(t* +T) > z(t*--nor)ws. 由 
式 (2) 可 得 
i(t) > z(t)(r —am4 —a1M) = eiz(t), t (n4 Dr, (11) 
a(t*) = (1 — p)z(t), t= (n Dr, 
ft [t*, t* + nar] 上 求解 (11)， 可 得 z(t* 十 n27) > maá(1— p)'2e?27, 由 此 可 知 
z( +T) 2 mall ~ petias > ma, 


TH. 令 f = jinf {x(t) > ma}， 于 是 当 t € (tt] HW. z(t) < mos Ba(t) 2 m. X Ttc 
(t*,t], E i(t) > oqa(t) 可 知 


z(t) 2 m4(1l — p)'atnsg(nana)air A ui: 


由 于 z 人 的 > ma， 因 此 对 于 t > E ERRIETA, “t> t, r(t) > mi. 

2) 如 不 是 脉冲 点 ， 则 t € [t,t*) Bb, z(t) > ma B.a(t*) = ma. fit € (n4 十 Dr， (ni 十 
L-- 1)r). Xft € (t*, (nj 十 1 十 1)7]，z(t) 的 取 值 有 两 种 可 能 

G) 对 所 有 的 t e (1*,(n4 十 1 十 1)7] 时 ， e(t) < ma， 则 类 似 于 情形 1) 可 以 证 明 一 定 存在 
RS € ((nj LU D)r, (nj - Ur 1)r +T] EBr) > ma. i= af {z(t) > ma}, 
则 te (t*, 0) tj, e(t) < ma; MŽ) =m. 4 Tte (t,0. 8 


z(t) > m4(1 E P)na+nae(na+ns+l)oar 全 Sm < m. 


HIT ali) > rn4a， 重 复 上 述 过 程 可 知 当 上 > 五 时，z(b > m. 
(说 ”存在 一 点 上 € (t*; (n 十 1 十 1)7]， 使 得 z(t) > m. 9£— inf {z(t) > ma}， 于 是 t € 
[t*, £) Ff, (t) < m4 Hel 2 ma. EAER 


z(t) > z(t*)e?10 > m4e^* > m. 


Hiz( > mm4， 重 复 上 述 过 程 可 知 ， 当 站 > 二 时 ，z(b) >m. 
综 上 所 述 ， 当 t > tı EF, z(t) > mı. 
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4 ”数值 模拟 


这 节 我 们 考虑 脉冲 周期 对 系统 的 影响 . Was 1, a = 2, a= 4, d=8, m=1,r=2, r= 
0.5, r2 = 0.8, p = 0.1, pı = 0.05, p; = 0.05， 数 值 模拟 如 图 1 至 图 3. 由 定理 2.1 和 定理 3.1 可 
知 ， 当 脉冲 周期 7 小 于 rmax 时 ， 害 虫 根除 周期 解 是 全 局 稳定 的 ， 取 由 = 01, uz = 0.2, T= 
0.2 如 图 1. 如 果 脉 冲 周期 7 超过 rmax， 由 定理 3.1 知 系统 (2) 是 持续 生存 的 ， 取 7 = 24 如 图 2. 
图 3 是 wi = 0.1, uz 从 0.001 增 加 到 8 时 的 分 支 图 ， 从 分 支 图 中 我 们 可 以 清楚 地 看 出 系统 (2) 随 
着 ua 的 增加 出 现 了 许多 丰富 的 动力 学 行为 ， 其 中 包括 : 拟 周期 震荡 、 倍 周期 分 支 、 对 称 破裂 
分 支 和 混沌 等 行为 ， 如 果 固 定 wa = 0.1, uj 从 0.001 增加 到 12 时 可 以 得 到 类 似 的 分 支 图 . 
































0.9 ^ -一 一 0.5 3 s 
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Br 时 间 序 列 图 ， 其 中 7 = 0.1 
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图 2: 时间 序 列 图 ， 其 中 7 = 24 
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图 3: 系统 (1) 和 系统 (2) 关于 参数 ua 的 分 支 图 ， 其 中 7 二 25 
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The Dynamics of a Predator-prey Model Based on 
Integrated Pest Management 


WANG Ling-shu'?, XU Rui?, FENG Guang-hui? 


(1- School of Mathematics and Statistics, Hebei University of 
Economics & Business, Shijiazhuang 050061; 2- Institute of Applied Mathematics, 
Mechanical Engineering College, Shijiazhuang 050003) 


Abstract: Based on the integrated pest management, a predator-prey system with impulsive effect is 
proposed and analyzed in this paper. By the Floquet theory and the comparison theorem of impulsive 
differential equation, it is proved that there exists a globally stable pest-eradication periodic solution 
when the impulsive period is less than some critical value. The sufficient condition for the permanence 
of the system is established. Numerical simulations illustrate the effect of impulsive on the dynamics 


of the system and the availability of our impulsive control strategy. 
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